Electronic doping of graphene by deposited transition metal atoms 
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We perform a phenomenological analysis of the problem of the electronic doping of a graphene 
sheet by deposited transition metal atoms, which aggregate in clusters. The sample is placed in a 
capacitor device such that the electronic doping of graphene can be varied by the application of a 
gate voltage and such that transport measurements can be performed via the application of a (much 
smaller) voltage along the graphene sample, as reported in the work of Pi et al. [Phys. Rev. B 
80, 075406 (2009)]. The analysis allows us to explain the thermodynamic properties of the device, 
such as the level of doping of graphene and the ionisation potential of the metal clusters in terms 
of the chemical interaction between graphene and the clusters. We are also able, by modelling the 
metallic clusters as perfect conducting spheres, to determine the scattering potential due to these 
clusters on the electronic carriers of graphene and hence the contribution of these clusters to the 
resistivity of the sample. The model presented is able to explain the measurements performed by 
Pi et al. on Pt-covered graphene samples at the lowest metallic coverages measured and we also 
present a theoretical argument based on the above model that explains why significant deviations 
from such a theory are observed at higher levels of coverage. 

PACS numbers: 73.22.Pr, 72.80.Vp, 73.30.+y 
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I. INTRODUCTION 

Graphene was discovered in late 2004^''^. This material 
is a one-atom thick sheet of carbon atoms, arranged in a 
honeycomb lattice. This structure is not a Bravais lattice 
and graphene is described in terms of a triangular lattice 
with a two-atom basis. A simple nearest-neighbour tight- 
binding approximation of the electronic Hamiltonian in 
graphene reveals that such lattice structure leads to a dis- 
persion relation that is linear around two specific points 
of the Brillouin zone. Since the Fermi level of graphene 
lies at these points, its quasi-particles behave in a contin- 
uum approximation as massless relativistic fermions with 
a speed of light equal to the Fermi- velocity « lO^ms"^ 
(see^^ and the recent review'*). 

The properties of graphene and its special geometry 
make it a very interesting candidate for applications in 
nano-electronics. Recent research has revealed other pos- 
sible applications, in solar cell technology^, in liquid crys- 
tal devices^, in single molecule sensors^, and in the fab- 
rication of nano-sized prototype transistors^. 

Transport measurements on graphene devices^ have 
become standard and can be performed under different 
doping conditions. Given the location of the Fermi level 
and the absence of a band gap between the valence and 
conduction bands in undoped graphene, one can contin- 
uously control the level of doping simply by the applica- 
tion of a gate voltage in a geometry where graphene acts 
as the upper (grounded) electrode of a capacitor. The 




FIG. 1: Capacitor device whose upper electrode is a single 
graphene sheet. The application of a gate voltage imposes a 
certain level of electron or hole doping on graphene, continu- 
ously increasing or lowering its Fermi level. The application 
of a potential difference between the upper contacts allows for 
transport measurements to be performed. 



lower electrode is composed of silicon, whereas the dielec- 
tric medium in between is Si02. Metal contacts placed 
on top of the graphene sheet allow for the realisation of 
transport measurements at different gate voltages, and 
hence at different levels of doping, with great flexibility. 
The system has an overall thickness of 6 « 300 nm (see 
figure [1]). 

The measurement of the transport properties using 
such devices can be used to determine the influence of 
different physical effects on both the AC and DC con- 
ductivities. One can investigate the influence of electron- 
electron interactions, of impurities, or of the presence 
of elastic ripples in the graphene sheet on the transport 
properties of this semi-metaP. 

It well known that undoped graphene, when analysed 
from the point of view of a self-consistent theory or the 
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renormalisation group, presents a finite conductivity with 
an universal value of 4e^/(7r/i)^"^^, regardless of the scat- 
tering mechanism that limits conductivity in graphene. 
The experimental measurements:-^'-- point to a some- 
what higher value for this quantity, equal to 4e^//i. This 
latter value is also obtained in studies of numerical diag- 
onalisation of graphene's tight-binding Hamiltonian with 
add-atoms acting as the source of disorde r^^i^^ . In the 
case of doped graphene, the behaviour of the conductivity 
markedly depends on the scattering mechanism that lim- 
its such quantity. It is therefore essential to clarify the na- 
ture of such a mechanism. The research community has 
held two opposing views, namely charged (Coulomb) or 
short-range scatterers^^, but recent experiments^ seem 
to show the latter mechanism as the prevailing one, even 
if it is agreed that charged scatterers also play a roleSS. 

In this paper, we will consider the contribution to the 
conductivity of one particular type of short-range disor- 
der, namely that induced by the deposition of transition 
metal (TM) atoms in graphene'^^. This type of disor- 
der is always present on devices such as those depicted 
in figure [1] due to the diffusion of metallic atoms from 
the contacts into the graphene sheet. The adsorption of 
graphene on TM surfaces has been extensively studied, 
both experimentally2i"— , as well as theoretically'^^"'^^. In 
TM surfaces for which the adsorption process (physisorp- 
tion) preserves the conical nature of the graphene bands 
close to the Dirac point (Al, Ag, Cu, An, Pt), the authors 
Qg5j36 j^ave shown that the levels of electron or hole dop- 
ing of graphene that they have found in their DFT studies 
can be explained by the relative value of the bulk work- 
functions of graphene and of that of the transition metal 
to which graphene is adsorbed. However, in order to ex- 
plain the electron-doping of graphene in cases where its 
work-function is lower than that of the transition metal 
(Ag, Cu), the authors invoked the existence of a chemical 
interaction between graphene and the underlying metal 
substrate, which plays a significant role in the formation 
of surface dipoles^l"— . The existence of such an interac- 
tion was confirmed in the experimental transport studies 
o^, performed on a graphene sheet where TM atoms 
were deposited, which was part of device such as that of 
figure [TJ The authors of this study have found that in the 
case of low coverage of graphene by Pt, the metal with 
the highest work function studied theoretically by^^i^, 
graphene is also electronically doped by Pt, becoming 
hole-doped at higher coverages. The authors stated that 
the high levels of electron-doping that they have found 
at low coverages were caused by an increased chemical 
interaction between graphene and the TM atoms, due 
to the short-distance (less than 3 A) between the two 
species (this distance is equal to 3.3 A in the full cov- 
erage regime). Furthermore, the AFM pictures obtained 
seem to show that the transition metal atoms aggregate 
in clusters at low coverage (see also^). It is neverthe- 
less unclear whether the proximity between the clusters 
and graphene is sufficient to justify the level of doping of 
graphene. 



The purpose of this paper is threefold. Firstly, we wish 
to introduce a framework that allow us to discuss the 
problem of charge doping of graphene by transition-metal 
clusters with generality from a thermodynamic point of 
view. This framework will be of a phenomenological na- 
ture and it will involve some simplifying assumptions, but 
it will already contain the main ingredients that will need 
to be considered in a more fundamental approach. Sec- 
ondly, the same type of phenomenological analysis will 
be extended to the problem of electronic scattering in 
graphene caused by the presence of the said clusters. We 
will show, following— that, despite the charged nature 
of the clusters, the scattering potential that they create 
is of a short-ranged nature (see also^^). The domain of 
validity of the semi-classical approximation of indepen- 
dent clusters that we are using is also discussed. Thirdly, 
these two elements of the theory will be used to interpret 
the above experiments from a quantitative point of view. 
This application of the theory will also serve to illustrate 
its overall limitations and we will provide physical ar- 
guments that show why more elaborate approaches are 
needed. 

The structure of this paper is as follows: in section 
im we will discuss the doping of graphene by metal clus- 
ters in a capacitor device based on a phenomenological 
model that treats each cluster as a perfectly metallic 
object kept at a constant potential dependent on the 
amount of charge in the cluster. The minimisation of 
the internal energy of the system at T = 0, subjected to 
overall charge conservation, will allows us to obtain the 
equilibrium conditions that determine the level of dop- 
ing of the graphene sheet. One can show, for equally 
charged clusters, that the level of doping can be writ- 
ten in terms of the bulk work-functions of the different 
components of the system, of the gate voltage applied to 
the device, and of a parameter that characterises the ef- 
fective chemical interaction between the graphene sheet 
and each individual cluster. The numerical value of this 
parameter is determined by two different contributions: 
the first contribution is due to the induced surface dipole 
of graphene and of the metallic clusters caused by the 
presence of the other components of the system; the sec- 
ond contribution is the correction to the Fermi energy 
of a cluster due to its finite size. Specialising to the 
case of spherical clusters, we can estimate the magni- 
tude of the chemical interaction using the measured val- 
ues by.?£ of the gate voltage that is necessary to apply to 
Pt and Ti-covered^ graphene samples to bring graphene 
to an uncharged state, where the conductivity is a min- 
imum. These values are dependent on the concentration 
of metallic atoms per unit cell as well as on which metallic 
element is deposited on graphene. In section IIIIl we will 
consider the form of the scattering potential created by a 
spherical cluster, following the model o^ and its contri- 
bution to the resistivity of the sample, within the First 
Born Approximation (FBA) and we will compare our re- 
sults with the transport measurements o^, performed 
on Pt-covered graphene samples at low coverage. We will 
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also show why the theory presented is not adequate to 
explain the measurements performed at higher coverages, 
both for Pt and Ti-covered samples. We will determine, 
using the above model, the linear dimension of the region 
in graphene where the charge donated by the cluster to 
this material is contained and show that, except for the 
lowest coverages considered in the experiments oS^, this 
quantity is comparable to the average distance between 
clusters, even for heavily-doped graphene. In section ITVl 
we will present our conclusions. In appendix \^ we will 
derive an expression for the cluster's ionisation potential 
that will be used in the main text, based on the same 
thermodynamic arguments that were used in section HIl 
Finally, in appendix |B] we will derive, using the method 
of images, the electrostatic contribution to the ionisation 
potential of a single spherical cluster, a result that will be 
shown to be in agreement with that of appendix [A) This 
derivation will also allow us to obtain the capacitance of 
the system composed of the metallic cluster and of the 
graphene plane, as well as the electrostatic potential due 
to a charged spherical cluster close to a grounded plane, 
a quantity that enters in the calculations performed in 
section IIIIl 



II. ELECTRONIC DOPING OF GRAPHENE BY 
DEPOSITED METAL CLUSTERS 

At T = 0, the internal energy of a composite system 
of k conductors can be written^i^, assuming that the 
electrons and the (immobile) ions of the different species 
interact with each other via the bare Coulomb interaction 
(i.e. one is including the contribution of the low lying 
electronic orbitals explicitly in the energy), as 



£ = 



G(iVi,.. 

4i: 



y, 47reo 



Pf i'^) P] 



(1) 



where the functional G{Ni, • • • , Nk) includes the kinetic, 
exchange and correlation energies of the electrons and 
where the second term includes the effect of the ion po- 
tential on the electrons and the Hartree energy of these 
electrons, with pf (r) = Pii'i') — pf"'' {f) being the plasma 
charge density. The indices i,j run over 1, . . . , fc. One can 
write pf (r) = p2{f) + ^Piif): where p"(r) is the charge 
density in the neutral ground state of each conductor 
and 5pi{r) is the excess charge density of that conductor 
due to charge exchange with the others. In particular, 
d'^rSpi{r) = Qi, the total unbalanced charged con- 
tained in conductor i. Using this decomposition, one can 
write the ground-state energy, up to a constant term, as 

£ = G(iVi,...,A^fe) + V / d'rVrir)Sp,{r) 

+ 1^ /" f ^d„, ^P^{r)5pJ{r') 



K Jv, 47reo \ r-r' 



(2) 



where V"(t") is the potential on conductor i due to itself 
and the other conductors, each in a neutral state. In 
a classical approximation, Spi{r) will be non-zero only 
close to the surface of the conductors, and we can write 
the above expression in a capacitor approximation: 

£ = G{Nl,...,Nk)+J2D^Q^ 

■i 

+ \T.^^'Q^Q^^ (3) 

where C^^^ is the inverse cross-capacitance between con- 
ductors i and j and Di = V"(r) is the surface dipole 
of conductor i^^, which is the average value of the elec- 
trostatic potential within that conductor. Note that the 
surface dipole of a conductor is computed with the re- 
maining conductors present, but in a neutral state. Thus, 
one expects that such surface dipoles will depend both on 
the geometry of each conductor and also on the presence 
of the other conductors if the distances between them are 
on the atomic scale. 

As stated above, the experiments of reference^^ were 
performed on a devices similar to that depicted in figure 
[1] with the transition metal atoms deposited by molecu- 
lar beam epitaxy (MBE) on the graphene sheet at dif- 
ferent coverages cs = Nam/nu of metallic atoms per 
unit cell of graphene, where Nam is the total number 
of deposited metallic atoms and Uu is the number of 
graphene's unit cells. In order to model such a device, 
we will assume that the atoms aggregate in identical 
clusters, which are randomly distributed above the area 
Ag — n-uAc of the graphene sheet, where Ac is the area of 
the graphene unit cell. We also assume that these clus- 
ters are all equally charged. The graphene sheet is kept 
at zero potential. At a distance b « 300 nm below it, one 
places a Si layer, with the space in between filled with 
Si02, a medium of permittivity e = 3.9 cq. The graphene 
sheet and the Si layer are connected to a battery such that 
a constant gate- voltage Vg is kept between them ('plane- 
capacitor model', see figure [5]). A single cluster-graphene 
subsystem is assumed to possess a joint capacitance Cs 
(which is computed for spherical clusters in appendix IB|) . 
The capacitance of the graphene/Si02/Si device is given 
by Csi = eAg/b. We assume that the cross-capacitance 
effects between different clusters are only due to the pres- 
ence of the grounded graphene plane, an assumption that 
is correct for small C5^^. With such assumptions, the in- 
ternal energy ^ can be written for this system, as 

£ = GiNs,...,Ns,Ng,Ns^)-eNcDsiNs-N|) 
- e Dg (Ng - A^O) - e Ds^ (A^5. - N^,) 

+ ^ (A^5 - + ^ {Ns. - Nl.f . (4) 

We have written the charge Qi of a given component i as 
Qi — — e{Ni — Ni), where Ns, Ng and Nsi are, respec- 
tively, the number of electrons in a cluster (all clusters 
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FIG. 2: Schematic representation of model that is considered 
in this paper. Transport measurements can be performed at 
different coverages as more atoms are deposited by MBE. 



Nam /name, where Hamc IS the number of metal atoms 
per cluster, which is equal to Uamc = zsVs/vs, with vs 
being the volume of the metallic unit cell and zs being 
the number of atoms in the unit cell (zg — 1 for Pt, 
= 2 for Ti). Expressing Nam hi terms of cs, which was 
introduced at the beginning of this section, one obtains 
for Nr 



Nr. 



(9) 



are equally charged), in graphene and in the Si layer in 
the equilibrium state in which these materials are in con- 
tact and exchange charge, and Ng, Ng and Ng^ are the 
same quantities in the uncharged state of these materi- 
als. Also, Ds, Dg and Dsi are the surface dipoles of each 
substance. 

The conditions of thermodynamic equilibrium are ob- 
tained through the minimisation of ^ subjected to the 
constraints that the overall charge of the system is zero 
and that the potential difference between the two elec- 
trodes of the battery is equal to Vg- The minimisation 
condition is thus given by 



-eVa dNs^ + /i {N, dNg + dNg + dNsi) , (5) 



since the charge transferred between the electrodes of the 
battery \& dq — —e dNsi, and where /i is the chemical po- 
tential of the system. One obtains from ([5]) the following 
equilibrium conditions 



= f^s + ^{Dg - Ds) 



e^Ns-Nl) 



Cs 



(6) 



= l^si + - Dsi) + eVa ■ 



where /ig — \Nj=Nsj with Nj being the number of 
charges in cluster j, Jig 



dG jj 
dN„ ' f^Si 



dG 
^Ns^ 



chemical potentials of the different components of the 
system in the absence of a dipole layer. These quantities 
are also called the internal contribution of a metal to 
its work-function^^. Note that ^ — Jig — eDg, i.e. the 
chemical potential of the system is equal to the (full) 
chemical potential of the subsystem kept at zero voltage. 

The two equilibrium conditions (|6l7p are not sufficient 
to determine the level of doping of the constituents of the 
system. These equations have to be supplemented with 
the neutrality condition for the overall system. One can 
write this condition as 



Nr ANs + AN„ + AN: 



Si 



0, 



(8) 



with ANs ^ Ns - N^, ANg ^ Ng ~ N^^ and ANs^ = 
Nsi — Ng^. The quantities that appear in ([5]) are related 
to the variations of the carrier density of the cluster, of 
graphene, and of the Si layer, by ANs = ^ns, ANg — 
Ag Sug and ANsi = Ag Snsi, where Vs is the volume of 
the cluster. The number of clusters is given by Nc = 



where we have expressed as the ratio between the area 
of the graphene sheet and the unit cell area. Substituting 
this formula in equation ([SJ and expressing ANs, ANg 
and ANsi in terms of the variations of the charge density 
of each media, one obtains 



csvs 
zsAc 



Sns + Sug + Snst 



0. 



(10) 



The equations (|6l7p , when written in terms of the vari- 
ations of density 6ns and Snsi, become 



Ms + <Dg - Ds) 



f^g — H'S^ "^y^g 

'Pa 



Cs 



Sns , 



(11) 



Jlsi + e{Dg ~ Dsi) -f eVc + — Snsr , (12) 



These two equations, which impose the equality of the 
so-called electro-chemical potentials^ between a metal- 
lic cluster and the graphene sheet, and between the 
graphene sheet and the Si layer, when supplemented by 
([To)) , are sufficient to determine the level of doping of 
graphene. However, we still need to relate the chemical 
potential Jig to 5ns, Jig to Sug and /J^i to 5nsi, in other 
words, we need the equation of state for the different com- 
ponents of the system. One writes Asp = Jig — Sp for 
the Fermi energy variation of the cluster, Ae^, = Jig — Sp 

Si _ 



are the — 



for the Fermi energy variation of graphene and Asp 
Psi ~ tl^^ Fermi energy variation of the Si layer, 

measured with respect to the uncharged ground state 
of each of these constituents. In the case of the clus- 
ters or of the Si layer, one has Sng ~ ps{£p) Aef, and 
Snsi ~ Psii^P) Aep since the density of states p{£f) is 
approximately constant for these materials at the Fermi 
level. Substituting these definitions in (fTTj) and (fT2|) and 
taking into account that the bulk work functions of the 
transition metal, of graphene, and of Si, are given by 
Wb = ei?s-£|, W^g = eD"g-el, and W^^ - eD^-ef^, 
one obtains 



Sns = 



eiADg~ADs)-Cs- 



A4 



1/Ps(4) 



e^Vs/Cs 



(13) 



Sn 



Si 



-W^ + e{ADs^-ADg-VG) 



Ae% 



-(14) 



where C,s — " the difference between the Fermi 
energy of the TM cluster and the Fermi energy of the 
bulk transition- metal and ADs = Ds — Db, ADg = 
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Dg — Dg and ADsi = Dsi — -0°, are the induced surface 
dipoles on each component of the system due to finite 
size effects and to the presence of the other components. 
One can estimate psii£p) ~ the result for a free 

two-dimensional electron gas, where m*g^ « me, is the 



electron's effective mass in Si. With b w 300 nm, one has 
Psii^p) ^ ^/i^'^b) ^iid one can neglect the first term in 
the denominator of (HH). One can thus write (113114^ as 



5ns 
5nsi 



Wb 
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A4 



l/ps{e%) + e^Vs/Cs 



(15) 
(16) 



where Ac = Cs + s{ADg — ADs) represents a cor- 
rection to the doping of the clusters due to their fi- 
nite size and to the induced surface dipoles, and Vq — 
(W^^ - W°)/e + AD Si - ADg. The quantity A^ can 
be interpreted as giving the overall magnitude of the ef- 
fective chemical interaction between the clusters and the 
graphene sheet. 

The density of states Pg{e) of graphene is zero at the 
Dirac point and one needs to consider its full functional 
form in that neighbourhood. It is approximately given 
by25 



(17) 



where t = 2.7 eV is the first-nearest neighbour hop- 
ping matrix element in graphene. Note that the pres- 
ence of impurities in graphene, either intrinsic or the 
deposited TM atoms themselves, will modify the den- 
sity of states given in (jl7p for high enough impurity 
concentrations^^. Integrating Pg{e) between the lower 

band limit Ep ~ \/ V^irt and the Fermi energy yields 
a result of two electrons per unitary cell. 

Integrating Pg{£) between Sp and pg, one obtains for 
dug the result 



Sui, ~ ± —= 



(A4f. 



(18) 



with the plus sign if Ae^ > and the minus sign other- 
wise. 

Substituting equations and ([H]) in (HUl), we 

finally obtain a second-degree equation for Asp 



± {Ae^pY + AAs^p-n^O., 



(19) 



with the plus sign if Asp > and negative sign otherwise, 
and where 



A 



V37rf2 

2 



csvs 



zs(l/ps(4)+e'Vs/Cs) e^b 
csvs {W° + A,-Wb) 
zs{llps{£p) + e^VslGs) 



+ '-^{Vg-Vo) 
eb 



(21) 



If we take the positive sign in equation (IT^ then a pos- 
itive solution exists if > 0. Conversely, if we take the 
negative sign in this equation, a negative solution exists 
if < 0. One can thus write for Ae^, the solution^ 



A4 = sign(r!) 




(22) 



One can see from equation (j22p that for a given con- 
centration cg, one can, through the application of a gate- 
voltage Vd such that = 0, bring the graphene sheet to 
its uncharged state, as Asp = 0. The gate voltage Vd 
can be determined from transport measurements on Pt 
or Ti-covered graphene^^, since the conductivity will dis- 
play the minimum characteristic of the Dirac point for 
that applied voltage. Likewise, the gate voltage Vq can 
be determined from the same measurements performed 
on the uncovered graphene samples, since it follows from 
equation ^ that for eg = 0, Ae^p = at Vb = Vb. 
Thus, the conductivity will also display the characteris- 
tic minimum at this applied voltage. Therefore, one can 
extract Ac from experiment. It is given by 



Ac = l^B - 



ezsA^jl/psjef) + e^Vs/CsKVp - Vp) 
e b csvs 



(23) 



The density of states at the Fermi level of bulk 
Pt or bulk Ti can be extracted from specific heat 
measurements^!^, through p{ep) = 37/(7rfcB)^, using 
the general result from Landau's Fermi liquid theory, 
where 7 is the linear coefficient for the dependence of 
the electronic specific heat cy on the temperature. Fi- 
nite size corrections to the bulk density of states may be 
estimated for a spherical cluster, using the results for a 
free-electron gas^, as psi^p) = pi^p) ~ -im* /{Snti^R), 
where m*p^ « 2m^, m'^^ w 3.15me^'^ is the electron's ef- 
fective mass in platinum or titanium and R is the radius 
of the cluster. Also, for a spherical cluster, Vs — 47ri?'^/3 
and Cs can be written as a power series on a parameter 
dependent on the cluster radius and on the distance L of 
its centre to the graphene sheet (see appendix [B|) . We 
take = 4.5 eV as ii^, Wl>^ = 5.64 eV and W^^ = 4.33 
eV for polycrystalline platinum and polycrystalline tita- 
nium, and Wsi = 4.6 eV^S,. We note that the model as 
defined contains two unknown parameters, namely R and 
L. 

Using these results, as well as the values of C5, Vq and 
Vd measured by^^ (raw data is a courtesy of Kawakami's 
group) and taking the radius of the cluster to be 0.6 nm 
for Pt^ and 0.188 nm for TiS and the distance from the 
centre of the cluster to the plane equal to 0.85 nm for Pt 
and 0.44 nm for Ti, one obtains for Ac the results given 
in tables HI [ll] and IIIII Note that one does not dispose 
of direct information (e.g. from AFM measurements) 
regarding the values of R and L. The values indicated 
above were chosen such as to provide agreement between 
the values of Ac and A^"* in tables HI |lT] and IIIII as well 
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CS 


Vd (V) 


R (nm) 


Ac (eV) 


AT (eV) 


PS 


0.025 


-11.4 


0.6 


2.47 


2.55 


-0.014 


0.071 


-29.0 


0.6 


2.49 


2.55 


-0.014 


0.127 


-46.0 


0.6 


2.37 


2.55 


-0.013 



TABLE I: Values of the concentration of Pt atoms per unit 
cell of graphene cs and applied voltages Vd corresponding to 
the minimum of conductivity for the Pt-1 sample studied in^^ 
(Vb = —1.94 V). The values of Ac were determined using for- 
mula (|23p . We have also computed the number of electrons ps 
per Pt atom at the Dirac point. This number was estimated 
in^'' as -0.014 e/Pt atom. 



cs 


Vd (V) 


R (nm) 


Ac (eV) 


A^* (eV) 


PS 


0.0065 


1.56 


0.6 


2.16 


2.55 


-0.011 


0.019 


-2.56 


0.6 


2.26 


2.55 


-0.012 


0.039 


-11.0 


0.6 


2.45 


2.55 


-0.014 


0.064 


-24.3 


0.6 


2.66 


2.55 


-0.016 



TABLE II: Values of the concentration of Pt atoms per unit 
cell of graphene cs and applied voltages Vd corresponding to 
the minimum of conductivity for the Pt-3 sample studied inS^ 
(Vo = 3.41 V). The values of A^ were determined using for- 
mula (|23p . We have also computed the number of electrons ps 
per Pt atom at the Dirac point. This number was estimated 
ir^ as -0.019 e/Pt atom. 



as between the theoretical and experimental asymptotic 
values for the contribution to the resistivity coming from 
the presence of the clusters^ (see section UlI)) . We have 
considered here and below the measurements made with 
samples Pt-1, Pt-3 and Ti-1 (in the notation of^^), since 
these samples, when uncovered, presented the smallest 
values of Vq measured, indicating a low level of intrinsic 
disorder. 

One can estimate the correction to the Fermi energy, 
due to the cluster's finite radius, from the free-electron 
gas result, as Cs ~ 8(m-)^^fl ' Pi'^p)- Using the re- 
sult quoted in the reference a^^i^^ for the induced dipole 
e {ADg - ADs) ~ 0.9 eV^?, one obtains for * the 
results presented in tables HI |TT] and IIIII for Pt and Ti, 
respectively. Thus, it is seen that the larger value of the 
chemical interaction Ac with respect to the case studied 
iii2^i2fi (particularly in platinum that has a larger DOS 
at the Fermi level) is due to a large shift of the Fermi 
energy of the clusters with respect to that of the bulk 
TM metal, caused by their finite radius. 

For a small cluster, the concept of work-function is 
ill-defined^°, as this quantity depends on the cluster's 
charge. One then speaks, respectively, of the cluster's 
ionisation potential if one is withdrawing an electron from 
a cluster at equilibrium, or of the cluster's electron affin- 
ity, if the electron is withdrawn from a negatively over- 
charged cluster. In appendix we compute the ioni- 
sation potential of a cluster based on a thermodynamic 



cs 


Vd (V) 


R (nm) 


Ac (eV) 


Af* (eV) 


ps 


0.0038 


-18.6 


0.19 


1.70 


1.83 


-0.179 


0.0077 


-41.0 


0.19 


1.89 


1.83 


-0.198 


0.0115 


-61.0 


0.19 


1.89 


1.83 


-0.198 


0.0153 


-73.5 


0.19 


1.71 


1.83 


-0.180 


0.0191 


-82.4 


0.19 


1.52 


1.83 


-0.161 



TABLE III: Values of the concentration of Ti atoms per unit 
cell of graphene cs and applied voltages Vd corresponding to 
the minimum of conductivity for the Tt-1 sample studied 
(Vo = —0.57 V). The values of Ac were determined using for- 
mula (|23p . We have also computed the number of electrons ps 
per Ti atom at the Dirac point. This number was estimated 
in^'^ as -0.174 e/Ti atom. 
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FIG. 3: Ionisation potential of cluster as a function of the 
applied gate voltage, for the Pt-1 sample, with coverages cs = 
0.025,0.071 and 0.127 ML (rgb). 

argument. We obtain from (jA4ll the result 

Is^W° + eADg + ^-Ae'p, (24) 

for the ionisation potential of the metallic cluster, where 
Asp is given by p2p . Substituting in (|24|) the parameters 
as computed in table U we have plotted in Figures [3l |4] 
and [5] the result ([M)) as function of the applied voltage, 
for the different coverages cs considered inM, for their 
Pt-1, Pt-3 and Ti-1 samples^-. In these plots, we have 
ignored the (unknown) constant eAD^. Nevertheless, 
such constant shift should be obtainable from a plot of 
the experimental ionisation potential, and so provide an 
estimate of eADg. 

It is also shown in appendix|X]that the electron affinity 
As of a metallic cluster is given by As = Is — e'^/Cs, 
with Wg = i(/s + As). Thus, As < V^g < Is- Since the 
transfer of an electron from the cluster to graphene would 
cost an energy Is — Wg > and, conversely, the transfer 
of an electron from graphene to the cluster would cost the 
same energy Wg — As > 0, one sees that the equilibrium 
state defined by equations ([TU)) to (IT^ is indeed a stable 
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FIG. 4: lonisation potential of cluster as a function of the 
applied gate voltage, for the Pt-3 sample, with coverages eg ~ 
0.0065, 0.019, 0.039 and 0.064 ML (rogb). 
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FIG. 5: lonisation potential of cluster as a function of the 
applied gate voltage, for the Ti-1 sample, with coverages eg — 
0.0038, 0.0077, 0.0115, 0.0153 and 0.0191 ML (romgb). 



one. The theory exposed in this section constitutes the 
main result of this paper. 



III. SCATTERING OF ELECTRONS BY THE 
METALLIC CLUSTERS AND ITS 
CONTRIBUTION TO THE RESISTIVITY IN 
THE FBA 

In the previous section, we have computed the level of 
doping of a graphene sheet due to the presence of metal- 
lic clusters. We have also computed the lonisation poten- 
tial of a single cluster. These properties are equilibrium 
properties. However, the experiments o:^ measured the 
dependence of the conductivity of graphene on the dop- 
ing induced by the metallic clusters and by the applied 
gate voltage. In order to describe such dependence, one 
needs to determine the scattering potential on individ- 
ual carriers due to the presence of the clusters. We will 
determine such a potential for spherical clusters, in an 
electrostatic approximation, which allows for the use of 



the method of images. Note that in such approximation, 
the graphene sheet is an equipotential surface. There- 
fore, the Coulomb interaction due to the surface charge 
distribution of the clusters is perfectly screened by the 
surface charge distribution that it is generated on the 
graphene sheet. However, such induced charge distribu- 
tion is spatially varying and it will thus correspond to 
a local variation of the Fermi level of graphene'*'*. Such 
variation will enter in the Dirac equation that describes 
the low-energy properties of graphene as a scattering po- 
tential and will give rise to a variation of the conductivity. 
Note that the local variation of the Fermi energy due to 
a cluster is not a parameter of the model as in22, but 
depends on the level of doping of graphene. 

The clusters, of radius R, are placed at a distance L > 
R above the graphene sheet. Each cluster possesses a 
charge Qs = —eANs = —^nR^eSns- We place the 
origin of the coordinate axis aligned with the centre of the 
spheres, such that the graphene sheet is located at z = 
—L. In the subspace z > —L, the electrostatic potential 
can be approximately described by the superposition 



V{r) = ^ i;(r - r-j) 



(25) 



of the potentials due to the individual clusters, located 
at r, = {xi,yi,0). 

In appendix iBl we will show how v{r) can be written 
in terms of a series of image charges, located in the clus- 
ter, and their images, located below the graphene plane. 
These charges depend on Qs, R and L through a re- 
cursion relation. The displacement field in the subspace 
z > —L is given by D = — eq Vy(r), whereas it is equal 
to Z? = —eSnsiBz in the space between the graphene 
sheet and the Si layer. The discontinuity of its normal 
component at z = — L determines the local density of 
charge ag[x,y) in the graphene sheet. In terms of the 
density of carriers Sng{x, y) = —cTg{x, y)/e, one has 

dng[x,y) = -onsi + - 



dz 



(26) 



where the derivative with respect to z is evaluated at the 
location of the graphene plane, z = —L. The spatial 
average of 6ng(x, y) is given by equation pop . Averaging 

Sns- 



dVjr) 
dz 



equation p6|l. we thus obtain 
Furthermore, one can write equation (|26p as 



ecsvs 



Sng{x,y) = Sug + 



dV{r) 



dz 



dV{r) 



dz 



-L. 



(27) 



The precise form of the local density of carriers depends 
on the location of the metallic clusters. However, assum- 
ing that one can treat clusters as independent entities, 
one has that in the neighbourhood of a given cluster, 
located at the origin of the coordinates, one can approx- 
imate dSll) by 



Sng{x,y) fa Sng + ^ 



dv[r) 



dz 



dv{r) 



dz 



(28) 
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where v{r) is given by (jB12p . Using the previous results, 
one can also estimate that 



dv{r) 



dz 



eazsAcN, 



Sns , 



(29) 



where A^c is the number of clusters. Since such a number 
is supposed to be very large, this term is negligible and 
one has that 



Sng{x,y) w Sug + 



cq dv{r) 



dz 



(30) 
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Substituting equations (|B12|) and (|B13I) for v{r) in ((50)) 
and expressing Qs in terms of Sns a-s above, one obtains 
the following series for dng{r)^ 



Sng{r) 



2 i?3 6ns - i?2 A"(l + A2") 
on„ ^ > — — 



3.g(A,l) 
1 



- (1-a: 



2n\2 



3/2 



(31) 



where A = (i — \/ L'^ — R"^) and where g(A, 1) is given 
by (IBTI) . 

The local variation of the carrier density 5ng{r) in 
graphene is related to the local variation of the Fermi 
energy Ae|,(r) through equation (|18p . Assuming that 
one is far from the neutrality point of graphene and that 
Ae|,(r) and Ae^, have the same sign, one has that the 
difference U{r) = Ae|, — Ae|,(r) is given approximately 
by 



U{r) 



TTt^ Ac 6ns 

2V35(A,1) I Aef, 
1 



_R2 x^'il + X^"") 
I ^ 

I n=l 



(1 - A2")2 



(L2 _ R2) 



3/2 



(32) 



where Asp is given by ([22]) and Sns is given by (flSl) . The 
function U{r) is the electron scattering potential due to a 
single cluster and depends, in this approximation, on the 
cluster carrier density Sns ^^nd also on the level of doping 
of graphene itself (through its dependence on | Asp |). 
Note that U{r) is attractive if Sns < (the cluster is 
doped with holes, as seen in experiment), as one would 
expect. The plot of U{r) is given, for different values of 
the gate voltage Vq, in figure [S] 

We can also compute from (PT|) . for later use, the size 
of the region in graphene that contains a charge —Qs 
of equal magnitude to that of an individual clusters, in 
the case in which graphene is electron-doped (i.e. for 
V > Vd)- We integrate equation (PT|) within the disk 
r < Rs, the region whose size we wish to calculate. We 
obtain after some cancellations, the following equation 



FIG. 6: Log-log plot of scattering potential due to clusters for 
Pt-1 sample with coverage cs = 0.025 and applied voltages 
Va = 0, 10, 20 V. The dashed line has slope —3 and is a guide 
to the eye. 



for Rs 



Rl 



4i^VL2 - J?2 I 
3g{X,l)Sng 



A"(1 + A2") 



. (1_A2»)2 



i?| +(L2-i?2) (l+A^ 



1/2 



(33) 



This quantity should be compared with the average dis- 
tance Rc between clusters, that can be simply defined 
through the relation ttRI = Ag/N^ = ^^^^^f^, i-e. in 
terms of the average area per cluster. If one defines the 
ratio ^ = Rs / Rc between these two quantities, one has, 
noting that the second fraction in the infinite sum of (j33p 
can be simply approximated by R'^^ 
is approximately given by 



if Rs > L, that ^ 



e = 




3/2 



A"(H-A2") 



4(A,z5i^)3/^9(A,l) 

1/3 



(1 - A2")2 




1/3 



(34) 



where the first term is a constant for fixed L and R. 
The validity of the independent cluster approximation, 
assumed above when passing from p7|) to (f28|) . depends 
on the condition ^ < 1— being fulfilled. 

The contribution of the metallic clusters to the resis- 
tivity of the sample of graphene can be easily computed 
within the FBA, once the scattering potential is known. 
Within the semi-classical theory'* based on the Boltz- 
mann equation, the contribution of the clusters to the 
conductivity of the sample is given by 

e^vphFTciikp) 



CTcl 



(35) 



where vp 
kp 



VV3Ae/2 



is the Fermi velocity in graphene, 
^-^f^ is the momentum of a quasi-particle at the 
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Fermi surface of doped graphene, measured with respect 
to the Dirac point, and Tci{kp) is the transport hfetime 
of a quasi-particle at the Fermi surface, due to the scat- 
tering with the metahic clusters. The expression above 
already accounts for the double spin and valley degener- 
acy (existence of two independent Dirac points). Equa- 
tion psp is known to apply as long as kpl,, 3> 1, where 
If, is the electron mean-free path. This condition holds 
in the diffusive regime where graphene is highly-doped 
(i.e. far away from the Dirac point) and for low-impurity 
concentration. The inverse of Tci{kp) can be computed 
in the FBA, by the application of Fermi's Golden-Rule 



1 



Tc/(fcf 



27r 



7V,5]|(fc'|[/(r)|fc)|2 
k' 

x{l-k -k) S{ek' - £k) , (36) 



where Nc is as above the number of clusters, i.e. the 
number of scattering centres, \ k \ = kp, e]^ = hvp \ k\, 

E]^' = hvp \ k' \, k, k are the unit vectors in the direction 
of k and k' and U{r) is given by (f32l) . 

In (Uni, one needs to take into account the spinorial 
nature of the wave- functions | A; ) , | fc' ) . The spinor 
Uk{f) — (flf^): which is normalized over the area Ag 
of the sample is given by 




(37) 



where tan^k = ky/kx- The ± signs stand for states with 
the same momentum and opposite energies relative to 
the Dirac point. The expression for u\^'{r) is entirely 
analogous. 

Substituting ([57)) and the analogous expression for 
Mk'(T') in (j36p . converting the summation over k' into an 
integral, performing the integral over k' using the delta 
function and expressing Nc/Ag = -^^^^fj--, we obtain 
the following result in terms of an angular integral over 
the scattering angle (, 



1 



Tcl{kp) 



'K'^csVsVp{dns)^R^ 
12zsA,kpgHX,l) 

^ A" 



2ir 



^ 1 

n=l 



X2n 



rrx^ 



(38) 

2 



The expression (p8| . as it stands, cannot be written in 
terms of elementary functions. However, the resulting 
integral is elementary if one can substitute the expo- 
nential functions in the infinite sum by 1, i.e. if their 
exponents are very small. The largest exponent is the 
one coming from the term with n = \ and is equal 
to 2fci?L sin((/)/2). Thus, this approximation is valid if 



kpL < 1 
condition as 



Since 5ng — ±- 



5ng |< 



one can write the above 
Taking L w 1 nm, one ob- 



tains from this condition that | Sug 10 — 10 e/cm . 
This condition is obeyed for all experimentally applied 




FIG. 7: Contribution to the conductivity of graphene sample 
due to Pt-clusters (Pt-1 sample) as function of the applied 
gate voltage, for the coverages cs = 0.025, 0.071 and 0.127 ML 
(rgb), following the prediction of equation (|39p . The values 
of the two fitting parameters, the cluster radius R = 0.6 nm 
and distance L — 0.85 nm between the cluster center and 
the graphene sheet, are as given in Table H] Also shown are 
the minimal value for the conductivity Ae'^/{nh) predicted 
by the SCBA and the value of 4e^//i measured b y^^i^° . Inset: 
contribution to the conductivity from the clusters as measured 
b3^. 



voltages. Thus, substituting the resulting expression for 
Tci{kp) in (p5|l . one obtains 



o-ciikp) 



6zsAc I fag I 
TT CS vs (Sns)"^ 



(39) 



where amin = 4e^/(7rft,) is the minimal conductivity of 
undoped graphene. Using the equations (P^ . ^TE\i and 
for Asp, Sug and 6ns in ([5^ . with the relevant pa- 
rameters in these equations taking the values as given in 
tables in we plot below (see figure [7]) the contribution to 
the conductivity of the graphene sample due to the clus- 
ters for the different coverages considered in2^ for their 
Pt-1 sample, as a function of the applied gate voltage Vg 
and as predicted by equation . Note the slight asym- 
metry of the curves with respect to the Dirac point, due 
to the variation of the carrier density Sns of the cluster 
with the applied gate voltage. 

One can also easily compute from ([5^ the contribution 
of the clusters to the mobility of the samples, defined as 
Mc; = o'ci/ie I Sug |). One obtains 



hires 2 



(40) 



where ps = vs5ns/zs is the number of electrons per TM 
atom or doping efficiency. 

The comparison of the result obtained in ([M)) with the 
experimental results of Pi et al. requires that we extract 
from their experimental data for the overall resistivity 
of the sample, the contribution coming solely from the 
metallic clusters, since there are other types of scatter- 
ers contributing to the resistivity even in an uncovered 
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sample, as discussed in the introduction. Therefore, we 
make the fohowing hypothesis regarding the dependence 
of the overall resistivity of a sample of graphene on the 
applied gate voltage Vg and on the metallic coverage cs 

p{VG,cs,ni) = Pimp{Sng,ni) + Cs PciiSug) 

+ PMsiSng,cs,ni) , (41) 

where rii is the concentration of intrinsic impurities rij in 
the sample. The function Pimp{Sng,ni) is the contribu- 
tion to the resistivity coming from the intrinsic impurities 
of the sample, which we take to be a sole function of the 
doping and of n^. This function can be extracted from 
the measurements done at zero coverage by expressing 
p{VG,0,ni) as a function of Sug, using equations (|18p 
and (|22|) with cs — 0- The second term is the contri- 
bution to the conductivity due to scattering by a single 
cluster and is therefore linear in cs- Equation (15^ is of 
this form, as we can always express 6ns in it in terms of 
Sug through (fT5|) and (fT8|) . However, this equation pre- 
dicts an infinite resistivity due to the clusters at the Dirac 
point, since it assumes the system to be in the diffusive 
regime, an assumption that fails close to the Dirac point, 
as discussed above (see alsc^'', where a similar situation 
occurs). Since the samples show a finite conductivity at 
the Dirac point and at finite coverage, we cannot take 
as it stands. Instead, we write for ad = p~i^, the 
ansatz 



Tcl 



QzsAc \6ng 
IT vs R'^ (Sns) 



(42) 



where ctJ?; is an extra contribution to the conductivity 
due to the clusters, which acts as an additional fitting 
parameter. Finally, PAfs('5?T.g, cg, n^) is the contribution 
to the resistivity due to multiple scattering events involv- 
ing the metallic clusters, be it multiple scattering by a 
single cluster, scattering events involving different clus- 
ters, or events involving clusters and intrinsic impurities 
in graphene, and is a general function of cs, Sug and n^. 

If one were to assume that pMsi^^g, cs, rii) were neg- 
ligible, the function ( p{Vg,cs, rii) - p(Vg, 0, rii) )/cs, ex- 
pressed in terms of Sug, would be independent of the 
coverage eg, i.e. it would be a universal curve. This 
is not the case for the metallic coverages considered by 
Pi et al, since these coverages are simply too large for 
multiple-scattering to be neglected, as will be shown be- 
low. In figure [SJ we perform a fitting of the theory to 
the results obtained with sample Pt-1 at the lowest cov- 
erage studied, but the objective of such fitting is merely 
to show that with the parameters characterising the clus- 
ters as given in table HI the theoretical and experimental 
results have the same order of magnitude and show the 
same asymptotic behaviour. The purple continuous plot 
represents the inverse of the function given by , with 
ctJ?; chosen so that the maximum of this curve and the 
maximum of the red-dashed curve coincide (this is the 
only free fitting parameter). We see that one is able to 
reproduce the asymptotic behaviour of the experimental 




200 300 



FIG. 8: Function [p{Vg, cs,ni)— p{VG,0,ni)]/cs , as measured 
for the Pt-1 sample, expressed in terms of the doping level 
Sug for the coverage eg = 0.025 (red-dashed curved), plotted 
against the inverse of (1421) (purple continuous curve). 
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FIG. 9: Ratio 5, as given by (|34|l . expressed in terms of the 
doping level Sug, for sample Pt-1, for the coverages cs = 
0.025, 0.071 and 0.127 ML (rgb). 



curve at large doping. The same asymptotic behaviour is 
also observed at higher coverages, but the curves deviate 
significantly from the universal curve hypothesis in the 
neighbourhood of the Dirac point. 

In order to understand the reason for the lack of agree- 
ment between the above theory and the experiments o^, 
we consider the behaviour of the ratio ^, introduced 
above, for the samples Pt-1 and Ti-1 (the behaviour ob- 
served for the sample Pt-3 is analogous to that of Pt-1). 
The plots are presented in figures and [TUl These plots 
indicate that, except for the lowest coverages, ^ > 1 in 
the whole range of doping displayed, and thus that the 
independent cluster approximation is unlikely to work 
for such high coverages. Note that this is purely a ge- 
ometric effect, caused by the small size of the clusters. 

One can also write Rr = 



= a 



Van 



where 



TTCS \ 27TCS 

a = 2.46 A is the length of the primitive cell of graphene. 
For the lowest coverage studied in the two Pt samples 
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likely to have a predictive power is below those coverages 
considered in the experiments of Pi et al. and thus one 
would need to repeat such experiments in these regimes 
in order to fully test such a theory. 
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FIG. 10: Ratio 5, as given by (|34|l . expressed in terms of 
the doping level Sug, for sample Ti-1, for the coverages eg — 
0.0038,0.0077,0.0115,0.0153 and 0.0191 ML (romgb). 



that we analysed, i.e. cs = 0.0065 ML for the Pt-3 sam- 
ple, with Uamc = 60 atoms, Rc ~ 12 nm. In the case 
of the Ti-1 sample, Uamc = 2 atoms and for the lowest 
coverage studied cs — 0.0038, i?c ~ 3 nm (for higher 
coverages, Rc is even smaller). Thus, in order to test 
the theory presented, the experiments performed would 
need to be repeated on samples presenting much lower 
concentrations of the deposited TM atoms (less than 1% 
for Pt-covered samples and less than 0.2% for Ti-covered 
samples). In addition, an appropriate characterisation 
of the clusters and their size distribution would also be 
required. With regard to the range of coverages stud- 
ied experimentally by^^, one should also note that if one 
takes the doping level of graphene to be Sug sa 10^^ cm^^, 
one has that kp « 10^ cm~^. Since Rc « 10 nm or less, 
kpRc ~ 1. This is yet another indication that a the- 
ory based on independent scattering centres is unlikely 
to work at this range of coverages. One should not ex- 
pect that the contribution of multiple scattering to the 
resistivity PMs{Sn-g,csTni) in (j4ip is a small quantity, in 
particular in the neighbourhood of the Dirac point^. 



IV. CONCLUSIONS 

In this paper, we performed a thermodynamic analysis 
of the problem of doping of graphene by TM clusters 
and computed the magnitude of the chemical interaction 
necessary to explain the electron doping of graphene by 
the transition metals Pt and Ti, the former having a 
bulk work function that is more than 1 eV larger than 
the work function of graphene. We have shown that 
the enhancement of such interaction with respect to 
the case studied in^^i^^ is due to the finite size of the 
TM clusters. We have also determined the scattering 
potential induced in a graphene sheet by spherical TM 
clusters and its contribution to the resistivity of the 
sample in the FBA. We have shown that regime of 
coverages for which the transport theory presented is 



Appendix A: Calculation of a cluster's ionisation 
potential based on a thermodynamic argument 

One can also determine the ionisation potential of a 
cluster from thermodynamic considerations. One starts 
by considering the expression for the energy (U) in the 
case of a single cluster {Nc = 1). Since the expression 
Q is valid for a system that is neutral, the extraction of 
a charge — e from the cluster requires this charge to be 
replaced in the graphene sheet. As discussed in appendix 
IbI the placement of such a charge in the graphene sheet 
can be pictured as the withdrawal of the image charge 
+e of the cluster charge — e, as the latter one is removed 
to infinity. Thus, one has from ^ that 

Is-Wg = £{Ns-l,Ng + l,Ns^)-£{Ns,Ng,Ns^) 
= -'Ps+'Pg + eiDs- Dg) 



(Al) 



Provided that the area of the graphene sheet is large, 
Wg ~ — /Ig + eDg , and thus 

2 2 

Is = -Ps + eDs-^ {Ns -Nl) + ^. (A2) 

Us ZUs 

One can consider the difference in the value of the ion- 
isation potential Is and the ionisation potential Ig in a 
situation where the cluster is placed very far away from 
the graphene sheet. One has 



Is 



-Ajls + e{Ds - D%) - 



1 

C~s 



1 



1 

C~s 



CI 



ANs 



(A3) 



Note that is not identical to Db, due to the finite 
size of the cluster. In appendix [B] we will show how the 
last two terms of (jA3p are obtained using the method of 
images to compute Is for the case of a spherical cluster. 

Finally, one can use the equilibrium condition © to 
write (jAip as 



IS^Wg 



2Cs 



eD„ 



2Cs 



(A4) 
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As stated above, one can explicitly compute the ca- 
pacitance of a spherical cluster of radius R, whose centre 
lies at a distance L > R, see appendix [BJ In this case, 
and in the limit R,L^ oo, with — R/2L finite, i.e. 
for graphene adsorbed on the bulk transition metal, as 
considered by^i^, Cs — > oo, and the transition metal 
ionisation potential and the work-function of graphene 
become equal at equilibrium, as one would expect (in this 
case, Is would just reduce to the TM work- function) . For 
a spherical metal cluster of finite radius, there is an extra 
contribution to its ionisation potential, coming from the 
effect of the image charge, already present in the case of 
an isolated spherical cluster—. This is the last term of 
the rhs of equation (|A4p . 

One can also compute the cluster's electron affinity As 
using the above argument. In this case, one is withdraw- 
ing an electron from an over-charged cluster and deliver- 
ing it to the graphene sheet. One has 

As~Wg - £{Ns,Ng,Ns^)-£{Ns + l,Ng~l,Ns^) 
= -'Ps+l^g + ^iDs- Dg) 



^(Ns-Nl) 



(A5) 



2Cs 2Cs 

where we have again used the equilibrium condition ^ . 
Note, in closing, that Wg = -^{Is + ^s), as one would 
expect. 



Appendix B: Calculation of the electrostatic 
contribution to the ionisation potential of a 
spherical cluster close to a grounded plane and of 
the electrostatic potential created by it 



displaced from z = d + R to z = oo, one needs to deter- 
mine the potential created by the presence of the sphere 
and of the plane at the position of the charge. Such po- 
tential can be determined by the method of images, as 
shown below. 

In the absence of a conducting plane, the solution of 
the problem is trivial. If qo is located at z ~ zq, one 
places an image charge of magnitude qi — ~qoR/ zq in 
the interior of the sphere at zi = R^/zq and an image 
charge q2 = Qs—qo—Qi at the centre of the sphere. These 
three charges guaranty that the surface of the sphere is 
equipotential and that the sphere has an overall charge 
Qs ~ Qo in it- In the presence of a grounded plane, these 
three charges no longer guaranty that the surface z — 
~L is an equipotential. We therefore take q2 to have 
an arbitrary value for the moment and consider three 
image charges (73 = —go at Z3 — — (2i + zq), (74 — —qi 
at Z4 = —{2L + zi) and (75 = —52 at Z5 = —2L, located 
below the plane. These three charges will guaranty that 
the plane z = —L is an equipotential. However, the 
surface of the sphere is no longer an equipotential. We 
therefore place three image charges in the interior of the 



sphere, qe = -(73 



NT 



— go 



2L+Z0 



Elt Zq — " I 



2L+Z0 ■■ 



2L+Z1 



^7 - - 



2L+Z1 



and 



-951 



at zs = = -g. In order 



97 = 
qs = 

to balance the potential at the surface of the plane, we 
now need to place three image charges below the plane, 
followed by three image charges inside the sphere and so 
on ad infinitum. The arguments above suggest that the 
recurrence relation between the charges inside the sphere 
is given by 



We will now compute the ionisation potential of a sys- 
tem composed by a single metallic spherical cluster and 
a grounded graphene plane using the method of images, 
by considering the electrostatic work necessary to extract 
a single electron from the cluster. These considerations 
will also allow us to write an explicit expression for the 
system's capacitance and for the electrostatic potential 
created by the cluster in the upper-half space, i.e. the 
quantity v{r), introduced in section Hill necessary to de- 
termine the scattering potential of carriers in graphene 
due to the presence of the spherical cluster. 

As above, the sphere contains a total charge Qs and 
we will extract a charge qo from it, leaving a total charge 
Qs ~ Qo in it- The charge go is the elementary charge 
— e. The ionisation potential of this system is the energy 
necessary to displace qo from a distance d away from 
the surface of the sphere to infinity. The distance d at 
which one begins to perform work to extract the charge 
qo is a regularisation parameter necessary to take into 
account the singular nature of the Coulomb interaction, 
but one can also interpret it physically as being the dis- 
tance beyond which quantum corrections to the Coulomb 
law become negligible. 

The conditions of the problem are as described in sec- 
tion mil In order to determine the force on qo as it is 



Q6(n+l)+a — QOn+a 2L+Z6„+a 
^6(n+l)+Q — 2L+Z6„ + c 



(Bl) 



where n > and a = 0,1,2. The recurrence relation 
between the charges located inside the sphere and be- 
low the plane is simpler, qsn+s+a = -Qon+a, zon+s+a = 

— {2L + ZQn+a)- 

The charge on the sphere is given by 



is = Qo + Ql + Q2 + Qo + Q7 + Qs + Q12 + Q13 + 914 



E 

n— 0,a 



(B2) 



This equation determines the charge (72 in terms of Qs 
and qo- Let us rewrite the recursion relation above in 
a shghtly different form. We define tt" ~ ^en+a, 
^Tl^ and 1] 



^ < i. We have 



'n+l 



W. 



n+l 



" l+w° 

nl_ 



(B3) 



We now define u" 



1 65 



we have = and 



It is easy to see from (jB3[) that 



Adding 



V 



n 1+10° 
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the two equations and using the recursion relation for w" , 
we have 



+ < 



II 
V 



(B4) 



This is a hnear recursion relation, with solution ~ 
A" s" + yl"s", where s± are the solutions of the 
quadratic equation ?7(sj_ + 1) — s± = 0, which are given by 

s± — zL ^-^2r) — ' ^ith s+s_ = 1 and s+ > 1,S- < 1. 
For definitiveness, we will call s_ = A, s+ = A^^. In this 
case, we have 



where we have used the recursion relation between the 
image charges on the plane and those on the sphere, and 
where the notation [n, a] indicates that the term n = 
0, a = (potential created by go) is absent from the first 
sum. The notation Vim{z, zq) indicates that the potential 
depends on zq through its dependence on the position 
and magnitude of the image charges. Expressing qen+a 
in terms of u" and z^n+a in terms of and using the 
rescaled variable w = z/2L, we have that Vim{w,WQ) is 
given by 



Vimiw,Wo) = 



1 



oo 



1 



= A1 A-" + A"_ A" 



(B5) 



where we have made use of the relation t](\ ^ + A) = 1. 
The values of the coefficients are determined from the 



oo 

E 



1 



<(u; + l + 



(BIO) 



Wo 



known values of Uq , Wq , namely Ug = ^ > 
zo/2L, ul = -Qo^^, "'o = £ and = 

One obtains Al = J- ^ — U^^(w,wq) ^ ^ 



Substituting the recursion relations given by (jB5l) above, 
we obtain for Uimiw^wo) = StteoL Vi,„(w, wq) 



1 



j^l _ 1 wo+ri\ 

+ QoV'^ A^i — A' — 



A WQ+riX 
qov A-i-A 



qoV A-i-A' 

and 



[n,a] 



yl^A~"(u; + r]X) + A^A"(u; + -qX^^) 



— — 7 , — ,_x ^ ■ Substituting these relations 

q2 A ^ — A' — 132 A ^ — A ^ 



and equation (IB5[) in equation (|B2[) . one obtains for q2 
the following result 

X (.g(A,0-r'5(A,r'))5"'(A,l),(B6) 



where ^ 



■i«o+r)A 



t«o+r(A 



^ < 1 andsa 



9i\0 



00 



^ 1 

ri=l 



(B7) 



^ yl^A-"(w + r;A-i) + A°LX"{w + rjX) ' 

(Bll) 

The potential v{r) created by the cluster on the sub- 
space z > —L, when its charge is equal to Qs, which was 
introduced in section Hill can also be computed in a man- 
ner analogous to Vim{z, zq). In this case one sets qo = 0, 
as in the calculation of the capacitance, in the recursion 
relation (IB5p . One is now interested in the dependence 
of v{r) both on z and on the radial coordinate r along 
the xy plane. 

This potential is given by 



Note that one can use this result to compute the ca- 
pacitance of the system composed by the sphere and the 
plane. For such a calculation, one takes qo = (or 
Wo Tj, i.e. ^ — >■ A). The charge of the sphere is Qs 
and the potential at its surface is given by Vs 



v{r) 



^ oo 



Q6n+2 



12 



'76n+2 



4;t£o J?, ' 



Hence, the capacitance Cs = Qs/Vs is, using (IB6[) with 
go = 0, given by 



[r2 + (z + 2L + Z6„+2)2]V2 



, (B12) 



Cs = 47reoi? (1 - A^) .g(A, 1)/A . 



(B8) 



where r — \/ + is the distance in the graphene sheet 
and where q%n+2 and zq„+2 are given by 



In terms of 77, the first few terms of this series are Cs 



Aire 



1 



rj- 



l-r;2 



(B13) 



The potential created by the image charges along the 
z-axis (z > R) is given by 



Vim{z, Zq) 



1 



E 



q^n+c 



47reo I '—'^ Z - Z^n+a 



E 



q^n+c 



Z -|- 2L -|- Zgn-l-o 



(B9) 



One can now use (IBllI) to compute the ionisation po- 
tential of the system cluster-graphene plane. This quan- 
tity is equal to the work of the external force necessary 
to transport the charge go from z = i?-|-c?uptoz = cxi 
quasi-statically, i.e. Is{d) = J^j^^ dzo Fext{zo), where 

P (-y \ — 10 9Ui^ I 

Substituting the values obtained above in (IB11[) 
and performing the derivative of lAim{w, wq) with respect 
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tow at w = wq. one obtains the rather lengthy expression 
for fe.t{wo) - ' ' 



/ext(wo) =qO^ 



{wo + ?7A-i)A-" - {wo + r]X)X" 



n=Q 



[ {wo + 77A-i)2A-" - {wo + ?7A)2A" ]2 



90 J2 

n=\ 
oo 

n=l 



(wo +?7A)A-" - {wo+7]X-^)X'- 



90 



E 



92 ^ 

ri=l 
oo 

92 ^ 

ri=l 



[ {wo + 


iixyx-^ - 

1 1 


- (wo + ??A-i)2A"12 


(t^o 


+ ryA-i)A- 


" - (wo + 77A)A" 


[(wo + 


r\X-^) {wo 


+ 77A) (A-« - A")]2 


{wo 


+ r7A)A-" 


- (wq +??A-i)A" 


[ {wo + 


r?A-i) (wo 


+ ?7A)(A-" - A")]2 




A-" 


-A" 


[{■Wo + 


7yA-i)A-" 


- (wo +ryA)A"]2 




A-" 


- A" 


[{wo + 


?7A)A-" - 


(wo + 77A-i)A«]2 ■ 



(B14) 



This expression has to be integrated so as 
to obtain the ionisation potential Is{d) = 
no^^ Q^dwofe.t{wo), where 7 = d/2L. It 
is more convenient to express this integral in terms of 
the variable ^ introduced above. The result is 



Is{d) = 



4 



1/5 



di{i-o 



+ 



1 X ^ /\ 

73 z^TY 



E 

00 



A" (1 - CA2") 

(l-^2;^2n)2 

A"(l-^^iA2") 



9gg(A,i) 

87reo\/l2~~R2 J^^^ 
9o 



47reoi? J 1/5 



1 ^ A 

'?2 



d?92(A,0 
A" (1 - A2") 



£2 ^ (1 _£-1^2 



^ A" (1 - A2") 
(B15) 



where 92(A,^) is given in equation (jB6[) . where we have 
also introduced (7(A, ^) and where 5 ■ 



> 1 and 



>)+7+'?A 

(5 < A~^, being equal to it in when d = 0. The first term 
in this expression can be integrated using the following 
identity 



1 



d 



2(1-0 

A"(l -^A2") 



{(l-0'[5(A,e')+r'.9(A,r')]}=(B16) 



00 



A"(l~C"iA2") 



^ (l-^2^2«)2 £3^ (l_£-2^2n) 



The second term involves a trivial integral. As for the 
third term, it can also be integrated if one notes that the 
following identity holds 



d 



{(l-0[.9(A,0-r\9(A,r')]}= (B17) 

2n\ 1 00 \ri ^^ _ \2n\ 



A" (1-^ 

71~TA2")2 £2 (1 _ £-1^2„)2 



^ A" (1 A2") l_ ^ X' (1 - A2") 
(1 _ £A2")2 ~ 7? 71 



Substituting these results and (jB6|) in (|B15[) . performing 
the resulting integrals and putting go = — e, one obtains 
after some trivial manipulations, the final result 

Is{d) 



lQTTeo{L + R + d) 



VL^ - i?2 



iTTeo[L + R + d-VL^-R^]^ 
X [g{X,6') +6-^g{X ,S-')] 
_ e2 VL2 - i?2 g(A, 1) 
~ 2T:eo [d^ + 2{R + L){R. + d)] 
{e^ + eQs)g-HX,l) 



47reo [ i + + - R2 ] 

X [g{X,5) -S-'g{X ,S-')] 

e2 y/L^-R'^ 5-'(A,l) 



47reo [ i + i? + d - /L2~7-r2 ]2 
X [g(A,(5)-5-i5(A,ri)]^ (B18) 

In the limit L — > 00, we are left with an isolated 
charged cluster. We obtain from (jB18|) the known re- 
sult 



R 



R 



47reo V 2d(d + 2i?) 
47reo(i? + d) ' 



R + d 2(i? + d)2 

(B19) 



The only term in (|BT8)) that is singular in the limit d — > 
is the first term of the series g{X,5'^) {6 — >■ A~^ in 
this limit). Moreover, for finite d, this term is equal to 
the corresponding singular term in (jB19[) . Thus, $(d) = 
Is{d) — I%{d) is a regular function in the limit d — > 0. It 
equals 

with Cg = AneoR and where Cs is given by (jBSp . Since 
Qs — —eANs, we see that the result (|B20[) is equal 
to the last two terms of (|A2p . which correspond to the 
electrostatic contribution to Is, the only one considered 
here. This justifies the statement made in appendix 1X1 
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uum theory such as the one presented above is applicable. 
The relative agreement obtained may be accidental, see^. 
The experimental results allow one to place a limit 
on the minimal distance between the clusters and the 
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L — R = 2.5 A, although the results for are only very 
weakly dependent on L — R, unlike what happens with R. 
The result of referencea^i^ for the induced dipoles was 
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In order to achieve faster numerical convergence^, 
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